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We studied the effect of magnetic doping with Gd atoms on the superconducting properties of 
amorphous Mo7oGe3o films. We observed that in uniform films deposited on amorphous Ge, the 
pair-breaking strength per impurity strongly decreases with film thickness initially and saturates at 
a finite value in films with thickness below the spin-orbit scattering length. The variation is likely 
caused by surface induced magnetic anisotropy and is consistent with the fermionic mechanism of 
superconductivity suppression. In thin films deposited on SiN the pair-breaking strength becomes 
zero. Possible reasons for this anomalous response are discussed. The morphological distinctions 
between the films of the two types were identified using atomic force microscopy with a carbon 
nanotube tip. 

PACS numbers: 74.48.Na, 74.25.Dw, 74.40.+k 



Understanding physical processes related to localized 
magnetic moments is particularly important for low- 
dimensional systems since such moments can form spon- 
taneously on surfaces and interfaces of nominally non- 
magnetic materials. The formation of localized magnetic 
moments is well known in semiconductor heterostruc- 
tures and devices, where they are carried by structural 
defects with unpaired electrons [1]. Localized magnetic 
moments were recently detected on the surface of a nor- 
mal metal |2]; in superconducting systems, they are be- 
lieved to be responsible for several unusual effects such 
as 1/f noise in SQUIDs and qubits Q and an anoma- 
lous magnetic field enhancement of a critical current 
in nanowires [4]. The origin of spontaneously formed 
magnetic moments often remains unknown; on the other 
hand, their effects can be probed by magnetic moments 
that are introduced intentionally. 

Here we study the effect of intentional magnetic dop- 
ing on transport properties of ultra-thin MoGe films that 
undergo a superconductor-insulator transition (SIT) 
[Tq| . The mechanism of the SIT remains an impor- 
tant unresolved problem in condensed matter physics. 
In general, there are several distinct physical processes 
that may lead to the SIT. Within the fermionic mecha- 
nism, Cooper pairing is locally suppressed by disorder- 
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enhanced electron-electron repulsion 
fermionic theories predict that the pair-breaking strength 
of magnetic impurities does not change with increasing 
disorder or decreasing film thickness 13|, LL4j . Experimen- 



tally, magnetic doping was studied in quench-condensed 
Pb[l5| and Pb-Bi films [l6|. In the latter case, behavior 
consistent with the fermionic theory was observed rela- 
tively far from the SIT. Several bosonic mechanisms were 
proposed for the critical regime of the SIT. In these mod- 
els, Cooper pairs are preserved across the transition but 
coherence in the films is lost due to vortex proliferation 
[TtJ, disorder-induced Cooper pair localization [Til fl9| , or 
fluctuations of the superfluid order parameter [20[ . While 



the models cited in Ref. [18l-l20j differ in their detail mi- 
croscopic mechanisms, they all predict the appearance of 
a spatially inhomogeneous superconducting state. The 
emergence of this state was observed in numerical sim- 
ulations [2l| and was recently detected experimentally 
[22| . Possible effects of magnetic pair-breaking within 
the bosonic models have not yet been analyzed theoreti- 
cally. 

The amorphous MoGe system is particular suitable for 
studying magnetic doping. This is the only known sys- 
tem where suppression of the critical temperature can be 
explained by the fermionic theory in all range of films 
thicknesses. Moreover, this can be achieved with the 
constrained theory, which assumes that effective electron- 
phonon coupling is not affected by disorder or film thick- 
ness. MoGe films with this property need to be deposited 
on a substrate covered with an underlayer of amorphous 
Ge that helps to maintain constant bulk resistivity of 
the film and ensures its homogeneity [23|. On the other 
hand, a missing Ge underlayer makes it possible to ob- 
tain and test an inhomogeneous superconducting state. 
We selected Gd as a magnetic dopant because its mag- 
netic moment is carried by a half-filled /-shell and does 
not depend on the host material. 

The critical temperature of amorphous Mo^Geioo-^ 
alloys depends on the particular value of x. In the first 
stage of our study, we used co-sputtering from three inde- 
pendently controlled guns with Mo, Ge, and Gd targets 
to fabricate a series of thick Mo^Geioo-^-Gd films with 
varying Gd content and x in the range 50-80 at. %. From 
transport measurements on these films we found that the 
alloy with x ~ 70 is the most suitable for the Gd dop- 
ing. In this alloy the superconductivity is completely 
suppressed when 6.5 at. % of Gd is added; at lower Gd 
content we detected a single-step superconducting tran- 
sition in R(T) curves. 

We fabricated two series of Mo7oGe3o films. Films of 
the A-series were deposited on a Si substrate covered with 
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FIG. 1: (A) Sheet resistance versus temperature for a series of 
amorphous Mo7oGe3o films deposited on SiN substrates (solid 
lines). The dashed lines indicate films doped with Gd; ar- 
rows indicate the correspondence between undoped and doped 
films. (B) Critical temperature versus sheet resistance for two 
series of undoped Mo7oGe3o films. The red solid line indicates 
a fit to the fermionic model. (C) Conductance of the insu- 
lating film 12 as a function of temperature on a logarithmic 
scale. The solid red line is a linear fit. 




a 60-nm thick layer of SiN grown by chemical vacuum 
deposition. For the B-series, prior to the deposition of 
MoGe film, a 3-nm thick underlayer of amorphous Ge 
was deposited. For oxidation protection the films of both 
series were covered by a 3-nm thick layer of Ge. 

In Fig. 1A, we show the temperature dependence 
of the sheet resistance for undoped Mo7oGe3o films de- 
posited on SiN (A-series). Within the studied tempera- 
ture range (down to 0.3 K) the system undergoes a di- 
rect SIT with no intermediate metallic phase. As shown 
in Fig. 1C, in the insulating regime, conductance has 
the logarithmic temperature correction arising due to the 
weak localization and electron-electron interaction con- 
tributions. Qualitatively similar suppression of supercon- 
ductivity was observed for films deposited on Ge under- 
layer. Figure IB shows the mean- field T c (defined at the 
middle of the transition) as a function of sheet resistance 
for the two series of films. As expected we found that 
the T c for the B-series can be well fitted by the fermionic 
model . However the T c of the A-series deviates from 
the model for the films with thickness below 1.5 nm. It 




FIG. 2: (A) 2D atomic force microscopy image of the surface 
of 1-nm thick MoGe film deposited on SiN (A- type). (B) 
2D AFM image of a similar film deposited on SiN covered 
with 3-nm thick Ge underlayer (B-type). Size of the images 
is 200 x 200 nm 2 



is interesting to note that the critical sheet resistance 
of the A-series (~ 5 kQ) is close to the universal sheet 
resistance R q = h/Ae 2 = 6.45 kTt predicted within the 
"dirty boson" model [ItJ The difference between A- and 
B-series cannot be explained by change in the dielectric 
constant of the substrate. SiN has lower dielectric con- 
stant than a-Ge; therefore, electron-electron interactions 
in this system have worse screening and suppression of 
T c would be expected at lower values of sheet resistance 
than in B-series. 

Looking for a possible structural effect of the Ge un- 
derlayer, we inspected surface morphology of several 
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FIG. 3: (A) 3D atomic force microscopy image of the sur- 
face of 1-nm thick MoGe film deposited on SiN (A- type). (B) 
3D AFM image of a similar film deposited on SiN covered 
with 3-nm thick Ge underlayer (B-type). (C) The difference 
A A between the ID Fourier profiles of the AFM images, nor- 
malized by the average profile, A, for the A-type and B-type 
films. 



test samples with an atomic force microscope (AFM) 
equipped with a carbon nanotube tip that provided 2 
nm lateral resolution. Typical 2D images are shown in 
Fig. 2A,B and 3D AFM images in Fig3A,B. To extract 
quantitative information about the lateral scale of the 
topographical features, a spatial Fourier transform was 
performed on the AFM data. To reduce noise, the result- 
ing data were averaged over different directions in Fourier 
space, yielding a ID profile of A vs. k for both the SiN 
and Ge underlayer samples. In Fig. 3C we show the dif- 
ference between the ID Fourier profiles of the A-type and 
B-type films shown in Fig. 3 A and Fig. 3B, normalized 
to the average of the A-type and B-type profiles. 

The overabundance of intermediate spatial frequencies 
for the type-A sample indicates that the Ge underlayer 
smoothes the surface and suppresses topographical fea- 
tures with a characteristic lateral scale of about 15 nm. 
This smoothing of the surface is also evident in com- 
paring the AFM images. On the other hand, the aver- 
age surface roughness, which characterizes the height of 
the topographical features, is 0.4-0.5 nm for both A-type 
and B-type films. Evidently, the missing Ge underlayer 



introduces inhomogeneities that are not strong enough 
to form a granular structure. This is also evident from 
transport measurements. We see from Fig. 1A that even 
for thinnest films the superconducting transition remains 
sharp with a well-defined T c . There is no tail in R(T) as 
it is typically observed in granular materials [25|. Nor- 
mal state properties of the films also do not indicate the 
presence of strong inhomogeneities. From the theory of 
weak localization (eq. 4.47a in Ref. (24|) we estimated 
that the dephasing length at T=0.3 K is about 80 nm 
for our least resistive insulating film. Since we do not see 
any sign of the insulating behavior down to T=0.3 K, 
the one-electron localization length in our films should 
be larger (probably much larger) than 80 nm and thus 
cover many random "hills" and "valleys" of the films' 
morphological profile. 

The inset to Fig. 4A shows an arrangement of sam- 
ples and targets used for fabrication of the MoGe films 
doped with Gd. The deposition was carried out by co- 
sputtering from two guns. Several substrates were po- 
sitioned approximately at the same distance from the 
composite MoGe target, but at varying distances from 
the Gd target. For each position, the deposition rate was 
calibrated by profile measurement of a test thick film; the 
thickness of films was controlled by the deposition time. 
Films within each series were fabricated in the same run 
under the same vacuum and deposition conditions. They 
have the same thickness but systematically varying Gd 
content. 

Temperature dependence of sheet resistance for several 
representative Gd-doped films is shown as dashed lines 
in Fig. 1A. The magnetic doping simply shifts the su- 
perconducting transition, leaving its width and normal 
state resistance essentially unchanged. Figure 4 displays 
the critical temperature versus Gd content. Errors in 
the Gd content originate from uncertainty in the time of 
the deposition, deposition rate and positioning of a sam- 
ple holder inside of the chamber. In addition, we found 
that thin A-type films with the same nominal thickness 
deposited in the same run revealed random variations 
in T c . It is interesting to note that this effect was not 
detected in the B-type films deposited on Ge. In A-type 
films random deviations of T c from the average value were 
always accompanied by the change in the normal state 
R s of a film; in fact, there was no uncertainty in T c vs 
R s relation. The error in T c resulting from this effect is 
indicated by vertical error bars in Fig. 4B. It was esti- 
mated from measurements on several series of undoped 
films deposited in the same conditions as the doped ones. 

In Fig. 4B we showed with open circulus the depen- 
dence of the T c on Gd content for the thick MoGe films 
fabricated by three-gun (Mo,Ge and Gd) deposition in 
the extended range of doping . The dependence can be 
fitted by the Abrikosov-Gor'kov (AG) theory [26]. The 
critical concentration of Gd is 6.5 at. %; the correspond- 
ing volume critical concentration is n c = 3 x 10 21 cm -3 . 



The rest of the data were obtained with the two-gun de- 
position. The AG theory predicts that at low doping, T c 
behaves as ^(T c q — T c ) = 7ra/4. The total pair-breaking 
strength, a, is related to the pair-breaking strength per 
impurity as a = a p n pi where n p is the concentration of 
impurities. A linear suppression of T c is expected, and 
indeed was observed experimentally. The parameter a p 
computed from the linear fit to the data, is plotted as a 
function of the film resistance in Fig. 5. 

We found that with decreasing film thickness the pair- 
breaking strength in MoGe films deposited on Ge drops 
by about three times initially and saturates in films with 
thickness below 1.5 nm. From the fermionic theories we 
expect that a p « const; however this conclusion is made 
under assumption that the exchange coupling between a 
localized spin and conduction electrons doesn't change 
with decreasing film thickness or increasing disorder. 

The behavior of a p in MoGe films appears to be qual- 
itatively similar to the reduction of the Kondo contribu- 
tion in thin films of normal metals doped with magnetic 
Extensive studies of this effect revealed that 
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FIG. 4: (A) Critical temperature versus Gd content for films 
with the Ge underlayer. The inset indicates the arrangements 
of samples and targets in the sputtering chamber for two- 
gun co-sputtering (B) Critical temperature versus Gd content 
for films with indicated thickness deposited on the SiN sub- 
strates. The dashed line is a fit of the data in the extended 
Gd percentage range to the Abrikosov-Gor'kov theory. 
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FIG. 5: The pair-breaking strength per impurity, a p , as a 
function of the film resistance for films deposited on SiN and 
amorphous Ge 



it is stronger for impurities with integer spin [IK and de- 
pends on surface roughness . The effect was explained 
in terms of spin-orbit induced magnetic anisotropy for 
magnetic impurities in proximity to the film surface [lo| . 
The theory predicts that close to the surface the effective 
spin of a magnetic impurity is reduced. 

Both the Kondo effect and magnetic pair-breaking de- 
pend on total impurity spin and exchange interaction 
between this spin and conduction electrons. From the 
known diffusion coefficient D = 0.5 cm 2 / s [23f| and spin- 
orbit scattering time r so = 5 x 10 -14 s [3jJ we can es- 
timate the average spin-orbit scattering length in MoGe 
as i so = nm. Experimentally, the t so coin- 

cides with the film thickness below which the saturation 
of a p takes place. This observation suggests that the 
pair-breaking strength of a Gd atom is reduced when it 
is located within i so from the surface of a film; the growth 
of a p in thicker films corresponds to the increasing frac- 
tion of Gd atoms with bulk-like surrounding. In other 
words, we have a gradual transition form anisotropic to 
isotropic exchange. 

For all MoGe films deposited on Ge we found that the 
suppression of the superconductivity by the magnetic im- 
purities and thickness reduction are additive processes. 
A magnetic impurity introduced into a sup erconductor 
suppresses the order parameter locally [32|. The local 
suppression of the order parameter with decreasing film 
thickness is also a feature of the fermionic mechanism of 
the SIT. In this regard, the additivity of the two pro- 
cesses even very close to the critical point of the SIT is 
consistent with the fermionic mechanism of the supercon- 
ductivity suppression. Moreover, the a p « const relation 
that we found in our thinnest films agrees with the spe- 
cific prediction made within the fermionic model. 

Let us now discuss how the magnetic doping affects 
thin A- type films. As shown in Fig. 4B, in the film with 
nominal thickness of 1 nm, the pair-breaking strength at 
low doping becomes zero; adding magnetic impurities to 
the film does not change its T c . One possibility for this 
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anomalous response is that the spin-orbit induced mag- 
netic anisotropy gets stronger in films deposited on SiN, 
because it has larger semiconductor gap. It is also possi- 
ble that the anomalous response to the magnetic doping 
is related to the enhanced roughness of the thin A-type 
films, which may result in inhomogeneous superconduct- 
ing state. Spacial non- uniformity of the order parameter 
is a common ingredient of the bosonic models. Analysis 
of the effect of magnetic impurities within these models 
can perhaps explain our findings. 

In summary, we have studied the effect of magnetic 
doping on superconducting Mo7oGe3o films. For uniform 
films deposited on amorphous Ge, the suppression of su- 
perconductivity is consistent with the fermionic mecha- 
nism. In thin films deposited on SiN, the pair-breaking 
strength becomes zero. Further analysis is needed to ex- 
plain this anomalous response. 
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